Heat conductivity in linear mixing systems 
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We present analytical and numerical results on the heat conduction in a linear mixing system. 
In particular we consider a quasi one dimensional channel with triangular scatterers with internal 
angles irrational multiples of vr and we show that the system obeys Fourier law of heat conduction. 
Therefore deterministic diffusion and normal heat transport which are usually associated to full 
hyperbolicity, actually take place in systems without exponential instability. 
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Given a particular classical, many-body Hamilto- 
nian system, neither phenomenological nor fundamental 
transport theory can predict whether or not this specific 
Hamiltonian system yields an energy transport governed 
by the Fourier heat law Q . Heat flow is universally pre- 
sumed to obey a simple diffusion equation which can be 
regarded as the continuum limit of a discrete random 
walk. In consequence, transport theory requires that the 
underlying deterministic dynamics yield a truly random 
process. Therefore, it is not mere idle curiosity to won- 
der what class, if any, of many-body systems satisfy the 
necessary stringent requirements. Moreover, it now be- 
comes increasingly meaningful to seek example of many- 
body systems which, using dynamics alone, can be shown 
to obey the Fourier heat law. A large number of pa- 
pers have recently approached this problem, mainly via 
numerical simulations [^-[Tl|]. Leaving aside, for the pur- 
pose of the present paper, systems which conserve the to- 
tal momentum, the general picture which emerges is that 
positive Lyapounov exponent is a sufficient condition to 
ensure Fourier heat law. In particular the paper |l0[ was 
precisely aimed at answering this question. Indeed in 
1^ the thermal conductivity was studied for a Lorentz 
channel- a quasi one dimensional billiard with circular 
scatterers - and it was shown to obey Fourier law. The 
conclusion at which the above numerical computations 
point out appears quite natural. Indeed modern ergodic 
theory tells us that for K-systems, a sequence of mea- 
surements with finite precision mimics a truly random 
sequence and therefore these systems appear precisely 
those deterministically random systems tacitly required 
by transport theory. On the other hand we do not have 
rigorous results and in spite of several efforts, the connec- 
tion between Lyapounov exponents, correlations decay 
and diffusive properties is still not completely clear. In a 
recent paper ||l2| a model has been presented which has 
zero Lyapounov exponent and yet it exhibits unbounded 
Gaussian diffusive behavior. Since diffusive behavior is at 



the root of normal heat transport then the above results 
constitutes a strong suggestion that normal heat conduc- 
tion can take place even without the strong requirement 
of exponential instability. If this would be the case, then 
it may turn out to be an important step in the general 
attempt to derive macroscopic statistical laws from the 
underlying deterministic dynamics. Indeed systems with 
zero Lyapounov exponent have zero algorithmic complex- 
ity and, at least in principle, are analytically solvable. 




FIG. 1. The geometry of the model. Particles move in the 
region outside the triangular scatterers. The x coordinate 
goes along the channel and y is perpendicular to it. The 
two heat reservoirs at temperatures Tl and Tn are indicated. 
The length of each cell is i = 3, the base of the triangles is 
a — 2.19 and the distance between the two parallel lines is 
d = 1.8. The geometry is then uniquely specified by assigning 
the internal angles 9 and (j) 

In this paper we consider a two dimensional billiard 
model which consists of two parallel lines of length L 
at distance d and a series of triangular scatterers (Fig. 
1). In this geometry, no particle can move between the 
two reservoirs without suffering elastic collisions with the 
triangles. Therefore this model is analogous to the one 
studied in |l^] with triangles instead of discs and the 
essential difference is that in the triangular model dis- 
cussed here the dynamical instability is linear and there- 
fore the Lyapounov exponent is zero. Strong numerical 
evidence has been recently given that the motion in- 
side a triangular billiard, with all angles irrational with 



1 



TT is mixing, without any time scale. Moreover, an area 
preserving map, which was derived as an approximation 
of the boundary map for the irrational triangle, when 
considered on the cylinder, shows a nice Gaussian dif- 
fusive behavior even though the Lyapounov exponent of 
the map is zero . It is therefore reasonable to expect 
that the motion inside the irrational polygonal area of 
Fig 1 is diffusive thus leading to normal conductivity. 

In the following we present careful numerical investiga- 
tions of the heat transport in the system of Fig 1 both by 
direct numerical simulation of energy flow for the system 
in contact with thermal baths as well as via Green-Kubo 
approach. Our results provide convincing evidence that, 
if the angles 9 and (j) are irrational multiples of tt, the 
system obeys Fourier law and the coefficients of thermal 
conductivity computed via the two approaches, coincide. 

In our computations we model the heat baths by 
stochastic kernels of Gaussian type, namely, the prob- 
ability distribution of velocities for particles coming out 
from the baths is 
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for Vx and Vy, respectively. 
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FIG. 2. Internal local temperature as a function of the 
rescaled cell number m/N for the irrational case with 
e = {^/2 - 1)tv/2 and (j) ^ 1. The total number TV of cells 
is: N ^10 (A), iV = 20 (□), N = AO (x), and iV = 80 (solid 
line). Here Tl ~ 1.1, Tr — 0.9, I — 3. Notice the quite good 
scaling behavior of the temperature field. 

The total length of the channel is L = Nl where N 
and / are the number and the length of the fundamental 
cells. For the irrational angles we take 6 = (\/2 — l)''''/2 
and (j) — 1. By increasing L, the number of particles per 
cell must be kept constant. However, since the particles 
do not interact one may consider the motion of a single 
particle over long times and then rescale the flux. 



We turn now to the definition of the two relevant quan- 
tities: the internal temperature and the heat flux. The 
temperature field at the stationary state is calculated fol- 
lowing the idea used in Ref |]lO|, namely, we divide the 
configuration space in slices {Ci}. The time spent by the 
particle within the slice in the jth visit is denoted by tj 
and the total number of crossings of a slice Ci during the 
simulation is M. The temperature is defined by 



Eft, ' 
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where Ej{Ci) is the kinetic energy at the jth crossing of 
the slice Ci. Since the energy changes only at collisions 
with the heat baths, we define the heat flux as 
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where {/S.E)k = — Eout is the change of energy at 
the fcth collision with the heat bath and Nc is the total 
number of such collisions which occur during time tc- 
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FIG. 3. Scaling behavior of the stationary heat flux J as 
a function of the system size for the irrational case of Fig. 1 
(•) and for the rational case (solid triangle) (see later in the 
text). A*' is the number of fundamental cells. The best-squares 
fltting gives a slope —0.99 ± 0.01 for the irrational case and 
—0.78 ± 0.01 for the rational one. 

For sufficiently long integration times, both the inter- 
nal temperature field and the heat flux reach a station- 
ary value. We have checked that the temperature profile 
obeys the law given in Ref | pO| . For small temperature 
differences AT, it is a linear function, as illustrated in 
Fig. 2: 
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In order to study the dependence of the heat flux on 
the system size, we need to consider larger and larger 
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systems, keeping the particles density constant. The cor- 
responding heat flux (for a density of one particle per 
unit length) is J = Lj, where j is the flux computed 
with a single particle simulation. In Fig. 3, we plot 
the heat flux J as a function of the system size N. For 
the irrational case, the best fit gives J = AN^^, with 
7 = 0.99 ± 0.01 and A = 0.015. The coefficient of ther- 
mal conductivity is therefore independent on N, which 
means that the Fourier law is obeyed and its numerical 
value is K = — — 0.225. 

The point under discussion here is very delicate and 
numerical experiments must be very accurate and reli- 
able in order to reach clear conclusions. We have there- 
fore checked the validity of our result by an independent 
approach, via a Green-Kubo type formalism, by study- 
ing the diffusive properties of our model isolated from 
thermal baths. 
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FIG. 4. (a) Initial temperature distribution I3{x) = 1/T{x) 
with /3(0) = l.l,/3(±L/2) = 0.9. (b)The rescaled heat flow 
^^'^^^^^ versus time t for different values of chain length for 
the irrational case. The dotted line has slope 0.225. 

In our numerical calculations, we follow Ref 1), 
namely, we consider the system with periodic boundary 
conditions and with an initial temperature distribution 
(3{x) given in Fig. 4(a). Then we calculate how the heat 
fiows from the half hotter part {L/4 < \x\ < L/2) of the 
system to the half colder part ( —L/4 < x < -l-L/4). At 
time t = we take a Maxwellian distribution of velocities, 
namely, P{v,^y) = exp(-w2 ,^/(2r(x)))/^27rT(x). If we 
denote by Q{t) the energy contained in the cold half part 
of the chain then, if the system obeys the Fourier law, 
the quantity {{AQ{t)) — {Q{t) — (9(0))) must increase 
linearly with time t, {AQ{t)) = {2KT^VP)t. Clearly the 
linear increase takes places only for times smaller than 



the sound transit time across channel. The numerical re- 
sults are shown in Fig. 4(b) where the dotted line, which 
fits the initial linear increase of the curve for N=40, has 
slope K = 0.225 thus indicating a very good agreement 
with simulations with thermal baths. 

Another important characteristic which is relevant for 
transport properties is the decay of the velocity auto- 
correlation function. In Fig 5 we show the decay of the 
absolute value of the normalized velocity autocorrelation 
C{t) — ^"'°^[,2^(Q°jj*''^ for the irrational case. This function 
decays as \C{t)\ - t~°' with a w 3/2. 




FIG. 5. The absolute value of the velocity autocorrelation 
function iC(t)j averaged over 6 x 10"" orbits initially with unit 
velocity amplitude but random directions. The solid line has 
slope —3/2. The dotted line shows the statistical errors. 

In Fig. 6 we show the diffusive behaviour of our model 
isolated from thermal baths. A clear diffusive behaviour 
(Ax^) = 2Dt for the irrational case is observed with dif- 
fusion coefficient D = 0.15. 

In fact, in our model, we can establish a connection 
between the thermal conductivity k and the diffusion co- 
efficient D. Indeed the average time a particle takes to 
travel from the left to the right bath and vice versa is: 
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then according to Eq.(|^), we can write the heat current 
as: 
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where P{v,T) = {At: / {2T:Tf/'^)v'^exp{-v'^ /2T) is the 
distribution for the modulus of velocity of particles com- 
ing out from the baths (see Eq. (1)). The thermal con- 
ductivity K — —Lj/VT is thus, 

n=^D. (7) 

We have numerically tested this formula. The thermal 
conductivity calculated from Figs. 2-3 is n = 0.225. On 
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the other hand the diffusion coefficient calculated from 
Fig. 6 is D = 0.15 which, according to the above rela- 
tion, gives K = 3D/2^ 0.225. 

As expected, a completely different behaviour is ob- 
tained when the angles and (j) are rational multiples of 
TT. The case with 9 = 7r/5 and 4> = ""/S, is shown in Fig. 
6 (triangles), and leads to a clear anomalous diffusive 
behaviour indicating the absence of Fourier law. Cor- 
respondingly, the simulations with the heat baths (solid 
triangles in Fig. 4) indicate a divergent behaviour of the 
coefficient of thermal conductivity k ~ N°^'^'^. 



to derive Fourier law. In this respect billiards like mod- 
els are very convenient since they are more suitable for 
analytical and numerical analysis. 
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FIG. 6. Diffusive properties of our model isolated from 
thermal baths. In the irrational case (Ax^) — O.SOSi^""^ 
(□); In the rational case, (Aa;^) = 0.082t^"® (A). 
Ax^ = {x{t) — a;(0))^. In the numerical calculation, 4 x 10^ 
particles are used. The particles are initially at j: = (in 
the center of the chain) and the initial velocities obey the 
Maxwell-Boltzman distribution at temperature T = 1. 

In conclusion, when all angles are irrational multiples 
of TT the model shown in Fig. 1 exhibits Fourier law 
of heat conduction together with nice diffusive proper- 
ties and the numerical value of the thermal conductivity 
computed via a Green-Kubo approach agrees with the 
one obtained by direct numerical simulations with ther- 
mal baths. However, when all angles are rational multiple 
of TT, the model shows abnormal diffusion and the heat 
conduction does not follow the Fourier law. 

One may argue that the model considered here is 
somehow artificial and far from realistic physical mod- 
els. However the problem discussed here is quite delicate 
and controversial and it is highly desirable to understand 
which dynamical properties are necessary and sufficient 



[1] 
it 
~ 

[4] 
[5] 
[6] 

[7] 



R. E. Peierls, Quantum theory of solids, Oxford Univer- 
sity Press, London (1955). 

F. Bonetto, J. L. Lebowitz and L. Ray Bellet 



phys/0002052 



math- 



G. Casati, J.Ford, F.Vivaldi and W.Visscher, Phys. Rev. 
Lett. 52, 1861 (1984). 

T. Prosen and M. Robnik, J.Phys.A: Math. Gen.25, 3449 

(1992) . 

H. Kaburaki and M. Machida, Phys. Lett. A 181, 85 

(1993) . 

S. Lepri, R. Livi and A. Politi, Phys. Rev. Lett. 78, 1896 
(1997); Euro phys. Lett, 43, 271 (1998); Physica D119, 



3541, (1999); |2ond-mat/0112193| ; 



[8] 
[9] 

[10] 
[11] 



B. Hu, B. Li and H. Zhao, Phys. Rev. E 57, 2992 (1998); 
A. Fillipov, B. Hu, B. Li, and A. Zehser, J. Phys. A: 
Math. Gen. 31, 7719 (1998). B. Hu, B. Li and H. Zhao, 
Phys.Rev. E 61, 3828(2000); B. Li, H. Zhao, and B. Hu, 
Phys. Rev. Lett. 86, 63 (2001); B. Li, H. Zhao, and B. 
Hu, Phys. Rev. Lett. 87 069402 (2001); B. Li, L. Wang, 
and B. Hu, Phys. Rev. Lett. 88, 23901 (2002). 

T. Hatano, Phys. Rev. E 59, Rl (1999); T. Prosen and 
D. K. CampbeU, Phys. Rev. Lett. 84, 2857 (2000). 

C. Giardina, R. Livi, A. Politi and M. Vassali, Phys. Rev. 
Lett. 84, 2144 (2000); O. V. Gendelman and A. V. Savin, 
Phys. Rev. Lett. 84, 2381 (2000). 

D. Alonso, R. Artuso, G. Casati, and I. Guarneri, Phys. 
Rev. Lett. 82 (1999) 1859. 

K. Aoki and D. Kusnezov, Phys. Rev. Lett. 86, 4029 
(2001); C. Mejia-Monasterio, H. Larralde, and F. 
Leyvraz, Phys. Rev. Lett. 86, 5417 (2001); P. L. Garrido, 
P. I. Hurtado, and B. Nadrowski, Phys. Rev. Lett. 86, 
5486 (2001); A V Savin, G P Tsironic, A V Zolotaryuk, 
Phys. Rev. Lett. 88, 154301 (2002); A. Dhar and D. 
Dhar, Phys. Rev. Lett. 82, 480 (1999 ); A. Dhar, ibid, 86 



[12] 
[13] 



1859 (2001); G Casati and T Prosen Jcond-mat/0203331 
G. Casati and T. Prosen, Phys.Rev. Lett. 85, 4261 
(2000). 

G. Casati, ad T. Prosen, Phys. Rev. Lett. 83, 4729 
(1999). 



4 



